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ABSTRACT: 

A scheduling method in a switch is disclosed, the switch including input ports (25.1 1-25.14) and 
output ports (25.5-25.54), each input port having one queue each for each said output port, 
thereby each input port having as many queues as a number of said output ports. Each input 
ports is represented as a node in a first set of vertices of a bipartite graph (23.1 ) and each said 
output ports is represented as a node in a second set of vertices of said bipartite graph (23.2). 
Each non-empty queue between any of said input ports and any of said output pods is 
represented as an edge between a vertex in said first set corresponding to said any of said input 
ports and a vertex in said second set corresponding to said any of said output ports, thereby 
creating a set of edges (23.3). An initial matching is performed by finding a subset of edges 
such that no vertex is incident to more than one edge in said subset (23.4). An unweighted 
shakeup matching for a hitherto unmatched vertex x in said first set to a vertex y in the second 
set is performed by performing an appropriate one of the following two steps, wherein if x 
matches with more than one of said vertices in the second set, vertex y being selected at 
random: adding new matching (x,y) if y is not part of any previous matching (23.9); replacing a 
matching involving y by a matching (x,y) if y is already part of a matching (z,y) wherein z is a 
vertex in said first set (23.1 1). A weighted shakeup method also considers the length of each 
queue. 
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(54) Shakeup technique to improve scheduling techniques for input-queued switches 



(57) A scheduling method in a switch is disclosed, 
the switch including input ports (25. 1 1 -25. 1 4) and output 
ports (25.5-25.54), each input port having one queue 
each for each said output port, thereby each input port 
having as many queues as a number of said output 
ports. Each input ports is represented as a node in a first 
set of vertices of a bipartite graph (23.1) and each said 
output ports is represented as a node in a second set of 
vertices of said bipartite graph (23.2). Each non-empty 
queue between any of said input ports and any of said 
output pods is represented as an edge between a vertex 
in said first set corresponding to said any of said input 
ports and a vertex in said second set corresponding to 
said any of said output ports, thereby creating a set of 
edges (23.3). An initial matching is performed by finding 
a subset of edges such that no vertex is incident to more 
than one edge in said subset (23.4). An unweighted 
shakeup matching for a hitherto unmatched vertex x in 
said first set to a vertex y in the second set is performed 
by performing an appropriate one of the following two 
steps, wherein if x matches with more than one of said 
vertices in the second set, vertex y being selected at 
random: adding new matching (x,y) if y is not part of any 
previous matching (23.9); replacing a matching involv- 
ing y by a matching (x,y) if y is already part of a matching 
(z,y) wherein z is a vertex in said first set (23.11). A 
weighted shakeup method also considers the length of 
each queue. 
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Description 

[0001] The present invention relates to a queuing system with multiple queues requiring multiple services, and more 
specifically to an improvement of a scheduling technique tar input-queued switches. Specifically, the present invention 
5 may be embodied in methods using an unweighted and weighted shakeup approach; method for queuing in a switch 
using an unweighted and a weighted shakeup approach; an unweighted shakeup switch and a weighted shakeup 
switch; and a network system using a combination of weighted and unweighted shakeup switches. 
[0002] For a better understanding of the present invention some graph-theoretic terminology is defined herein. 

10 • A bipartite graph, as shown in FIG. 1 , is an undirected graph for which the vertices can be partitioned into two sets 
1 .1-1 .4 and 1 .5-1 .8 such that an edge 1 .9 is located only between the two sets. A set of queues requiring a set of 
services can be represented by such a bipartite graph. The existence of an edge between two vertices indicates 
that an entry in a queue represented by the first vertex requires a service represented by the second vertex. 

• A matching is a subset of edges such that no vertex is incident to more than one of the edges in the subset. Finding 
is a matching corresponds to the selection of a partial permutation between queues and services requested. 

• A maximum (size) matching is a matching with the largest possible number of edges, n and m denote the number 
of vertices and edges respectively in an underlying graph representing the present state. For an N x N system 
representing N 2 queues requiring N services n = 2N and m < A£. 

• A maximum weight matching is a matching M that maximizes the sum of the weights of the edges in M. 

20 

[0003] An important example of such a queuing system is a router switch used in applications like the Internet. The 
present application uses a switch as an important example of a queuing system that can be represented using a bipartite 
graph. It should be noted that the present invention should not be limited to a switch and can be applied to any queuing 
system that can be represented using a bipartite graph. 

25 [0004] The growth of the Internet places continuously increasing demands on routers and switches. Switch designs 
capable of handling aggregate bandwidths of tens or hundreds of gigabits per second, with commensurate increase 
in the number of input and output ports are presently available. Examples of such switches are described in N. MeKe- 
own, M. Izzard, A. Mekkittikul, W. Ellersick, and M. Horowitz, "The Tiny Tera: A packet core switch' in Hot Interconnects 
IV, (Stanford University), pp. 161 -173, August 1996; and C. Partridge, P. P. Carvey, E. Burgess, I. Castineyra, T. Clarke, 

30 L Graham, M. Hathaway, R Herman, A. King, S. Kohalmi, T. Ma, J. Mcallen, T Mendez, W. C. Miliiken, R. Pettyjohn, 
J. Rokosz, J. Seeger, M. Sollins, S. Storch, B. Tober, G. D. Troxel, D. Waitzman, and S. Winterble. "A50-Gb/s IP router, 
' IEEE/ACM Transactions on Networking, vol. 6, no. 3, pp. 237-248, June 1998. 

[0005] In addition, such queues increasingly face the demands of supporting flows with requested quality-of-service 
(QoS), which may include bandwidth, latency, and jitter restrictions. Of course, all services must be supported in an 

35 efficient way, such that the switch will be able to operate at nearly 100% capacity. 

[0006] To support these increasing demands, a switch architecture that provides scalable bandwidth as the number 
of switch ports grows is promising. The examples described in the background section as well as the preferred em- 
bodiments described in the present specification assumes that the switching fabric is a crossbar switch that has N 
input ports and N output ports, N being any natural number. It should be noted that the present invention is not limited 

40 to such a switch and can be used for any switch including a switch with M input ports and N output ports, M and N 
being natural numbers. 

[0007] A crossbar switch is capable of routing an arbitrary permutation of fixed-sized packets (cells) in one time step. 
When Karol et al. illustrated the limitations of using a single queue per input port several research groups began to 
consider the use of N queues per switch-a technique known as virtual output queuing (VOQ). For more details on such 

45 limitations, See M J. Karol, M. G. Hluchyj, and S. P. Morgan, "Input versus output queuing on a space-division packet 
switch," IEEE Transactions on Communications, vol. COM-35, no. 12, pp. 1347-1356, December 1987. For more back- 
ground information, see also M. Karol, K. Bag, and H. Obara, "Improving the performance of input-queued ATM packet- 
switching," in Proceedings IEEE INEOCOM'92, pp. 110-115, 1992 and N. W. McKeown, Scheduling Algorithms for 
Input-Queued Cell Switches. PhD thesis, University of California at Berkeley, 1995. 

50 [0008] Intuitively, VOQ allows the switch to have greater flexibility when selecting a permutation. An example of a 
VOQ switch where N = 4 is shown in FIG. 2. For such VOQ switches, selecting appropriate permutations quickly is a 
fundamental problem. But what is the most appropriate permutation? One reasonable approach would be to route the 
most cells possible in the current time step. 

[0009] The above-mentioned problem is identical to calculating a Maximum Bipartite Matching (MBM), a graph prob- 
55 iem that has been studied in some detail (see Section IB.1(a)). However, the scheduling calculations involved in the 
implementation of conventional methods lead to severe bottleneck. This has led to a number of conventional heuristics 
that approximate the MBM. These heuristics include Anderson et al.'s PIM, McKeown et al.'s SLIP, and a number of 
other approaches that greedily find a maximal matching. See T E. Anderson, S. S. Owicki, J. B. Saxe, and C. P. Thacker, 
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■High-speed switch scheduling for local-area networks," ACM Transactions on Computer Systems, Vol. 11, no.4, pp. 
319-352, November 1993 and N. McKeown, R Varaiya, and J. Wairand, "Scheduling cells in an input-queued switch, 
■ Electronic Letters, vol. 29, no. 25, pp. 2174-2175, December 1993. 

[0010] A further drawback of MBM demonstrated by McKeown, is that it can be unstable even for traffic flows for 
5 which no input or output port is expected to handle a throughput greater than its maximum capacity. A queuing method 
is considered unstable when it can result in unbounded queue lengths. Further information on these drawbacks can 
be found in N. W. McKeown, Scheduling Algorithms for Input-Queued Cell Switches. PhD thesis, University of California 
at Berkeley, 1995. 

[0011] A flow is said to be admissible when no input or output port switch is loaded to its full capacity. It was later 
10 shown independently by Tassiulas and Ephremides; and McKeown, Anantharam, and Wairand that a generalization 
of the MBM approach, called Maximum Weight Bipartite Matching (MWBM), would be stable for ail admissible flows. 
See N. McKeown, V. Anantharam, and I. Wairand, "Achieving 100% throughput in an input-queued switch," in Pro- 
ceedings IEEE INFOCOM'96, (San Francisco, CA), pp. 296-302, March 1996. As with MBM, however, implementing 
full MWBM is problematic. This has resulted in a number of heuristics to approximate a MWBM, for example ILPF See 
15 A. Mekkittikul and N. McKeown, "A practical scheduling algorithm to achieve 1 00% throughput in input-queued switch- 
es" in Proceedings IEEE INFOCOM'98, (San Francisco, CA), pp. 792-799, March 1998. 

[0012] MBM and MWBM approaches may be suitable in different contexts. For best-effort service, where QoS is not 
considered, the MBM approach is quite reasonable. If there is no call-admission control for flows, then the ability to 
support all admissible flows will not be a particularly useful feature. In such a case, maximizing instantaneous bandwidth 
20 is a more reasonable design goal. On the other hand, if QoS is required and call-admission control is used, an approach 
like MWBM that is stable for ali admissible traffic is highly desirable. 

[001 3] While the techniques of the present invention choose a matching based only on the current state of the input 
queues, some related work considers the scheduling problem under slightly different assumptions. For example, Kam 
and Siu and Li and Ansari assume the scheduler has knowledge of flow requirements from the call-admission controller 

25 (CAC). See A. C. Kam and K. Y, Siu, "Linear complexity algorithms for bandwidth reservations and delay guarantees 
in input-queued switches with no speedup," in Proceedings of the 6th International Conference on Network Protocols, 
pp. 2-11, 1998 and S. Li and N. Ansari, "Input-queued switching with QoS guarantees," in Proceedings IEEE INFO- 
COM'99, vol. 3, pp. 1 1 52-1 1 59, March 1 999. Other conventional approaches calculate a matching while leveraging off 
the matching from the previous time step, for example, see A. C. Kam and K. Y, Siu, "Linear complexity algorithms for 

30 bandwidth reservations and delay guarantees in input-queued switches with no speedup," in Proceedings of the 6th 
International Conference on Network Protocols, pp. 2-11,1 998 and see L. Tassiulas, "Linear complexity algorithms for 
maximum throughput in radio networks and input queued switches," in Proceedings IEEE INFOCOM 98, vol. 2, pp. 
533-539, 1 998. The rationale in these cases is that the state of the input queues can only change by a small amount 
each time step; intuitively, the matching for time step i will very likely be a reasonable matching for time step i + 1 as well. 

35 [001 4] For a better understanding all scheduling approaches are partitioned into two classes; according to whether 
they approximate MBM (unweighted) or MWBM (weighted). These two classes are described in Sections IB.1 and IB. 
2, respectively. 

1B.1 Unweighted Approaches 

40 

[0015] The unweighted approaches to the scheduling problem do not take into account the size of the different 
queues, only whether they are empty or not. The switch can be modeled as a bipartite graph, where each input port 
is represented by a node in one set, each output port is represented by a node in the other set, and the existence of 
an edge between two vertices indicates that a cell exists to be routed between the corresponding input and output 

45 ports in the switch. Finding a matching corresponds to the selection of a permutation to route. In the present specification 
n and m are used to denote the number of vertices and edges respectively in an underlying graph representing the 
state of the switch. For an N x N switch, n = 2N and m =< A£. The representation of the switch as the bipartite graph 
defined above allows for succinct method descriptions. Almost all conventional scheduling heuristics aim to approxi- 
mate a maximum or maximum weight (see next section) matching. A related concept is that of a maximal matching, 

so to which it is not possible to add an edge while still maintaining a matching. The various conventional heuristic ap- 
proaches are described hereunder. 

IB.1 (a) MBM 

55 [001 6] A maximum bipartite matching approach is the first example of an unweighted approach. Computing an actual 
maximum size matching is now considered. The fastest method to compute a maximum size matching runs in 0(n 1/2 m) 
time worst-case. See E. A. Dinitz, "Algorithm for solution of a problem of maximum flow in networks with power esti- 
mation," Sower Math. Dokl, vol. 11 , pp. 1 277-1 280,1 970 and S. Even and R. E. Tarjan, "Network flow and testing graph 
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connectivity, ' SIAM Journal on Computing, vol. 4. pp. 507-518, 1975. 

[0017] The MBM approach has several drawbacks. It can lead to instability with an admissible but non-uniform traffic 
even on a 3 x 3 switch. See N. McKeown, V. Anantharam, and I. Walrand, "Achieving 100% throughput in an input- 
queued switch," in Proceedings IEEE INFOCOM'96, (San Francisco, CA), pp. 296-302, March 1 996. Another problem 
s with MBM is that the relevant methods are sequential in nature, requiring a centralized controller with knowledge of 
the state of each of the 0(n 2 ) queues at each time step. 

[0018] Because of these problems, several heuristics have been proposed that settle for a bipartite matching that is 
not maximum, but that can be calculated quickly (often in hardware) on real systems. See S E. Anderson, S. S. Owicki, 
J. B. Saxe, and C. P. Thacker, 'High-speed switch scheduling for local-area networks, 1 ACM Transactions on Computer 

10 Systems, Vol. 11, no.4, pp. 319-352, November 1993; M. AN and H. Nguyen, 'A neural network implementation of an 
input access scheme in a high-speed packet switch, 1 in Proc. of GLOBECOM 1989, pp. 1192-1196, 1989; T. X. Brown 
and K. H. Liu, "Neural network design of a Banyan network controller,' IEEE Journal on Selected Areas of Communi- 
cation, vol.8, no.8, pp. 1428-1438. October 1990; N. McKeown, P. Varaiya, and J. Walrand, "Scheduling cells in an 
input-queued switch," Electronic Letters, vol. 29, no. 25, pp. 2174-2175, December 1993 and T. R Troudet and S. M. 

15 Walters, "Neural network architecture for crossbar switch control," IEEE Transactions on Circuits and Systems, vol. 
38, no. 1, pp. 42-56, January 1991. 

IB. Kb) PIM 

20 [0019] Parallel Iterative Matching (PIM) is a heuristic, which will converge to a maximal matching in 0(log n) iterations 
on average. See T E. Anderson, S. S. Owicki, J. B. Saxe, and C. P. Thacker, "High-speed switch scheduling for local- 
area networks," ACM Transactions on Computer Systems, Vol. 1 1 , no. 4, pp. 31 9-352, November 1 993. At every iter- 
ation unmatched inputs propose to every output for which they have a queued cell. If an unmatched output receives 
any requests, it grants to one by selecting uniformly at random over all requests. If an input receives more than one 

25 grant it randomly selects one among them. 

IB. 1(c) RND 

[0020] The RND heuristic is similar to FIM. The difference is that every unmatched input /make only a single request 
30 to a randomly chosen output port / for which /" has a queued cell. 

IB. 1(d) RR-RND 

[0021] In the round-robin method RR-RND, scheduling for a given time slot t is determined by a fixed priority schedule. 
35 The input if with the highest priority chooses randomly among the outputs for which it has a cell queued. The input i 2 
with the next highest priority chooses next among the not-taken outputs, and so on. At the next time slot t + 1 the input- 
port priorities will be incremented by 1 , modulo N. The advantage of this priority scheme is that it can be easily imple- 
mented in hardware with a pipeline mechanism. 

40 IB. 1(e) SLIP 

[0022] SLIP was developed by McKeown and analyzed extensively in N. W. McKeown, "Scheduling Algorithms for 
Input-Queued Cell Switches", PhD thesis, University of California at Berkeley, 1 995. In SLIP, each input port and output 
port has a distinct priority wheel. An input sends a request to every output for which it has a queued cell. An output j 
45 chooses to accept a request from port i t based on its priority wheel; if the grant is accepted by /, the priority wheel at / 
will be set to (/+ 1 ) mod N. Input ports accept grants based on priority wheels, and advance their wheels in a similar 
fashion. 

[0023] The PIM, RND and SLIP heuristics fall into the class of iterative heuristics. Iterative heuristics are not guar- 
anteed to find a maximal matching in a single execution (iteration). Therefore repeating them more than once, may 
50 increase the size of the permutation. Of course there is a tradeoff with speed, since the switch fabric may have to 
operate at a higher speed if more than one iterations are required for each permutation. The number of iterations is a 
parameter of importance as shown subsequently (infra Section C). 

IB. 2 Weighted Approaches 

ss 

[0024] In the weighted approaches, the sizes of the queues affect the execution of the scheduling method. In the 
maximum weight bipartite matching problem the edges of the graph are assigned weights and the objective is to com- 
pute a matching of maximum weight. Maximum weight bipartite matching has the advantage that it has been proven 
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theoretically so be stable under admissible traffic where the traffic flow from input / to output j can be modeled as a 
sequence of independent and identically distributed random variables. See N. McKeown, V Anantharam, and I. Wal- 
rand, "Achieving 100% throughput in an input-queued switch," in Proceedings IEEE INFOCOM'96, (San Francisco, 
CA), pp. 296-302, March 1996 and A. Mekkfttikul and N. McKeown, 'A practical scheduling method to achieve 100% 

5 throughput in input-queued switches," in Proceedings IEEE INFOCOM'98, (San Francisco, CA), pp. 792-799, March 
1998. McKeown conjectures that MWBM is stable under all admissible flows, a conjecture that is also supported by 
the experiments. See N. W. McKeown, "Scheduling Algorithms for Input-Queued Cell Switches" PhD thesis, University 
of California at Berkeley, 1995. On the-other hand, MWBM methods are deemed impractical for hardware implemen- 
tation. For input-output pair (/, j) let qjj denote the occupancy of queue Q(/, j), i.e. the number of cells queued at input 

10 /with destination / 

IB.2(a) MWBM-LQF and MWBM-LPF 

[0025] Two distinct flavors of maximum weight matching schedulers have been studied, each defining the weights 

is on the edges of the bipartite graph in a different way. 

[0026] In the MWBM-LQF approach the weight Wg of an edge is set equal to c^. See N. McKeown, V. Anantharam, 
and I. Walrand, "Achieving 100% throughput in an input11 queued switch," in Proceedings IEEE INFOCOM'96, (San 
Francisco, CA), pp. 296-302, March 1996. In the MWBM-LPF method, first weights for the ports are defined. See A. 
Mekkittikul and N. McKeown, "A practical scheduling method to achieve 100% throughput in input-queued switches," 

20 in Proceedings IEEE INFOCOM'98, (San Francisco, CA), pp. 792-799, March 1998. In particular set Sj = X output^ 
for an input /. Analogously for an output /, Rj = X inputjqy. For the weighted bipartite matching instance define Wjj to be 
Sj + Rj if qjj > 0, and 0 otherwise. Both MWBM-LQF and MWBM-LW are stable under admissible traffic with independent 
arrival processes. See N. McKeown, V Anantharam, and I. Walrand, 'Achieving 100% throughput in an input-queued 
switch," in Proceedings IEEE INFO COM'96 t (San Francisco, CA), pp. 296-302, March 1996 and A. Mekkittikul and 

25 N. McKeown, "A practical scheduling algorithm to achieve 100% throughput in input-queued switches," in Proceedings 
IEEE INFOCOM'98, (San Francisco, CA), pp. 792-799, March 1998. The fastest method to compute a maximum-weight 
bipartite matching runs in 0(n 2 log n + nm) time worst-case. See R. G. Busaker and P. J. Gowen, "A procedure for 
determining minimal-cost flow network patterns," Tech. Rep. ORO-15, Operational Research Office, Johns Hopkins 
University, Baltimore MD, 1961; W. S. Jewell, "Optimal flow through networks," Tech. Rep. 8, Operations Research 

30 Center, MIT, Cambridge, MA, 1958; and M. L. Fredman and R. B. Tarjan, "Fibonacci heaps and their uses in improved 
network optimization algorithms," Journal of the ACM, vol. 34, pp. 596-615, 1 987. If the maximum weight of an edge 
is bounded by C, the cost-scaling method of Gabow and Tarjan runs in 0(n 1/2 mlog(nC)) time worst case. See also J. 
B. Orlin and R. K. Ahuja, "New scaling algorithms for the assignment and minimum mean cycle problems," Mathematical 
Programming, vol. 54, pp. 41 -56, 1992 and A. V. Goldberg and R. Kennedy, "Global price updates help," SIAM Journal 

35 on Discrete Mathematics, vol. 1 0, pp. 551 -572, 1 997 for alternative cost-scaling methods. 

IB.2(b) ILPF 

[0027] ILPF is a simple heuristic proposed in to approximate the solution found by MWBM-LPF. See A. Mekkittikul 
40 and N. McKeown, "A practical scheduling algorithm to achieve 100% throughput in input-queued switches," in Pro- 
ceedings IEEE INFOCOM'98, (San Francisco, GA), pp. 792-799, March 1998. The heuristic exploits the special as- 
signment of edge weights as a function of port weights. The inputs and outputs are sorted according to their S and R 
values. The sorted list of outputs is traversed in decreasing order and output /is matched to the input / with the largest 
Sj value such that (i) i is still unmatched and (ii) qg > 0. No justification is offered in why this heuristic may approximate 
45 the maximum weight matching found by MWBM-LPF. See A. Mekkittikul and N. McKeown, "A practical scheduling 
algorithm to achieve 100% throughput in input-queued switches," in Proceedings IEEE INFOCOM'98, (San Francisco, 
CA), pp. 792-799, March 1998. However, it is a well-known fact in matching theory that ILPF will find the optimum 
solution if the underlying bipartite graph is complete, i.e. an edge exists between all input-output pairs. See R. K. Ahuja, 
T. L. Magnanti, and J. B. Orlin, "Network Flows: Theory, Algorithms and Applications" Prentice Hall, 1993. 

so 

IB.2(c) RR-LOF 

[0028] This scheduler is similar to the RR-RND method described in Section IB.1(d). The difference is that each 
input /selects among the not taken outputs j the one with the maximum q^ value. 
55 [0029] It is an objective of the present invention to improve scheduling in high-throughput input-queued switches. 
The present invention uses a Shakeup technique, which can be used in conjunction with a number of the existing MBM 
and MWBM heuristics to improve their performance dramatically. The Shakeup technique assumes an initial matching 
is provided, and unmatched nodes have the ability to "knock out" existing matches. Intuitively, this provides a mecha- 
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nism to jump out of local maxima in the matching heuristic. 

[0030] To meet the objectives of the present invention, there is provided a computer-implemented unweighted queue- 
ing method for a system with M groups of N input queues each of said N input queues requesting a unique one of N 
services, wherein M and N are natural numbers, said method comprising representing each of said M groups as a 

s node in a first set of vertices of a bipartite graph; representing each of said N services as a node in a second set of 
vertices of said bipartite graph; representing each non-empty queue as an edge between a vertex in said first set 
corresponding to one of said M groups to which said each non-empty queue belongs to and a vertex in said second 
set corresponding to a service requested by said each non-empty queue; performing an initial matching by finding a 
subset of edges such that no vertex is incident to more than one edge in such subset; performing an unweighted 

10 shakeup matching for a hitherto unmatched vertex x in said first set to a vertex y in the second set by performing an 
appropriate one of following two steps, wherein if x matches with more than one of said vertices in the second set, 
vertex y being selected at random: adding new matching (x, y) if y is not part of any previous matching; replacing a 
matching involving y by a matching (x, y) if y is already part of a matching (z, y), wherein z is a vertex in said first set; 
and repeating for each unmatched vertex in said first set. 

is [0031] Another aspect of the present invention is a computer-implemented unweighted queuing method for a system 
with M group of N input queues each of said N input queues requesting a unique one of N services, wherein M and N 
are natural numbers, said method comprising: representing each of said M groups as a node in a first set of vertices 
of a bipartite graph; representing each of said N services as a node in a second set of vertices of said bipartite graph; 
representing each non-empty queue as an edge between a vertex in said first set corresponding to one of said M 

20 groups to which said each non-empty queue belongs to and a vertex in said second set corresponding to a service 
requested by said each non-empty queue; performing an initial matching by finding a subset of edges such that no 
vertex is incident to more than one edge in such subset, wherein if a vertex can be matched using more than one 
edges, one edge is selected from said more than one edges with a probability relative to the size of queues associated 
with each of said one or more edges; performing a weighted shakeup matching by performing an appropriate one of 

25 following two steps for a hitherto unmatched vertex x in said first set to a vertex y in the second set: if y is not part of 
any previous matching, said vertex y being selected with a probability: 



Jb=l 

35 wherein q^ is a length of a queue between vertices x and y; 

replacing a matching involving y by a matching (x,y) if y is already part of a matching (z,y), said replacement being 
done with a probability: 

40 q *y 

wherein z is a vertex in said first set; repeating for each unmatched vertex in said first set; and repeating for each 
unmatched vertex in said first set. 

45 [0032] Yet another aspect of the present invention is a method of unweighted routing in a switch comprising input 
ports and output ports, having cells to route between said input ports and said output ports, each input port having one 
queue each for each said output port, thereby each input port having as many queues as a number of said output 
ports, said method comprising: representing each said input ports as a node in a first set of vertices of a bipartite graph; 
representing each said output ports as a node in a second set of vertices of said bipartite graph; representing each 

50 non-empty queue between any of said input ports and any of said output ports as an edge between a vertex in said 
first set corresponding to said any of said input ports and a vertex in said second set corresponding to said any of said 
output ports, thereby creating a set of edges; performing an initial matching by finding a subset of edges such that no 
vertex is incident to more than one edge in said subset; performing an unweighted shakeup matching for a hitherto 
unmatched vertex x in said first set to a vertex y in the second set by performing an appropriate one of the following 

55 two steps, wherein if x matches with more than one of said vertices in the second set, vertex y being selected at random: 
adding new matching (x.y) rf y is not part of any previous matching; replacing a matching involving y by a matching (x, 
y) if y is already part of a matching (z,y) wherein z is a vertex in said first set; and repeating for each unmatched vertex 
in said first set. 
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[0033] Yet another aspect of the present invention is a method of weighted routing in a switch comprising input ports 
and output ports, having cells to route between said input ports and said output ports, each input port having one queue 
each for each said output port, thereby each input port having as many queues as a number of said output ports, said 
method comprising: representing each said input ports as a node in a first set of vertices of a bipartite graph; repre- 

s senting each said output ports as a node in a second set of vertices of said bipartite graph; representing each non- 
empty queue between any of said input ports and any of said output ports as an edge between a vertex in said first 
set corresponding to said any of said input ports and a vertex in said second set corresponding to said any of said 
output ports, thereby creating a set of edges; performing an initial matching by finding a subset of edges such that no 
vertex is incident to more than one edge in such subset, wherein if a vertex can be matched using more than one 

10 edges, one edge is selected from said more than one edges with a probability relative to a size of a queue associated 
with each of said one or more edges; performing a weighted shakeup matching by performing an appropriate one of 
following two steps for a hitherto unmatched vertex x in said first set to a vertex y in the second set: if y is not part of 
any previous matching, said vertex y being selected with a probability: 

75 

gxy 
N 

20 

wherein q^ is a length of a queue between vertices x and y; replacing a matching involving y by a matching (x,y) if y 
is already part of a matching (z,y), said replacement being done with a probability: 



wherein z is a vertex in said first set; repeating e for each unmatched vertex in said first set; repeating for each un- 
30 matched vertex in said first set. 

[0034] Preferably unweighted shakeup is performed for every alternate matching. 
[0035] Preferably weighted shakeup is performed for every alternate matching. 

[0036] Preferably unweighted shakeup is performed only if no increase in the size of matching is possible without 
shakeup. 

35 [0037] Preferably weighted shakeup is performed only if no increase in the size of matching is possible without 
shakeup. 

[0038] Yet another aspect of the present invention is an unweighted shakeup switch for use in a network, said switch 
comprising: input ports; output ports, wherein cells are routed between said input ports and said output ports; a graph 
generator that generates a bipartite graph with each said input ports as a node in a first set of vertices of the bipartite 

40 graph and each said output ports is a node in a second set of vertices of said bipartite graph and a cell to be routed 
between an input port and an output port is an edge between a vertex in said first set said second set; a matcher for 
performing an initial matching by finding a subset of edges such that no vertex is incident to more than one edge in 
such subset; an unweighted shakeup matcher for performing an unweighted shakeup matching for a hitherto un- 
matched vertex x in said first set to a vertex y in the second set, wherein said unweighted shakeup matcher adds new 

45 matching (x,y) if y is not part of any previous matching and replaces g a matching involving y by a matching (x,y) if y 
is already part of a matching (z,y). 

[0039] Still another aspect of the present invention is a weighted shakeup switch for use in a network, said switch 
comprising: input ports; output ports, wherein cells are routed between said input ports and said output ports; a graph 
generator that generates a bipartite graph with each said input ports as a node in a first set of vertices of the bipartite 

so graph and each said output ports is a node in a second set of vertices of said bipartite graph and a cell to be routed 
between an input port and an output port is an edge between a vertex in said first set said second set; a matcher for 
performing an initial matching by finding a subset of edges such that no vertex is incident to more than one edge in 
such subset; an weighted shakeup matcher for performing an weighted shakeup matching for a hitherto unmatched 
vertex x in said first set to a vertex y in the second set, wherein said unweighted shakeup matcher adds (x,y) if y is not 

55 part of any previous matching, said vertex y being selected with a probability: 
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N 

wherein is a length of a queue between vertices x and y, and replaces a matching involving y by a matching (x.y) 
if y is already part of a matching (z,y), said replacement being done with a probability: 

10 

is [0040] Preferably in the unweighted shakeup switch unweighted shakeup is performed for every alternate matching. 
[0041] Preferably in the weighted shakeup switch weighted shakeup is performed for every alternate matching. 
[0042] Preferably in the unweighted shakeup switch unweighted shakeup is performed only if no increase in the size 
of matching is possible without shakeup. 

[0043] Preferably in the weighted shakeup switch weighted shakeup is performed only if no increase in the size of 
20 matching is possible without shakeup. 

[0044] Still another aspect of the present invention is a network system comprising a 15 plurality of switches for 
routing messages, a first subset of the plurality of said switches being unweighted shakeup switches, a second subset 
of the plurality of said switches being weighted shakeup switches. 

[0045] The above objectives and advantages of the present invention will become more apparent by describing in 
25 detail preferred embodiments thereof with reference to the attached drawings in which: 

Fig. 1 is a bipartite graph abstraction of the scheduling problem in a switch; 
Fig. 2 is a crossbar switch that supports VOQ with N=4; 
Fig. 3 shows stability levels for PIM in regular and shake varieties; 
30 Fig. 4 shows stability levels for SLIP in regular, shake, and hybrid varieties; 

Fig. 5 shows stability levels for RND in regular, shake and hybrid varieties; 
Fig. 6 shows stability levels for RR-RND in regular and shake varieties; 
Fig. 7 shows stability levels for SHAKE; 

Fig. 8 shows stability levels for ILPF in regular, shake, and weighted shake varieties; 
35 Fig.9 shows stability levels for RR-LQF in regular, shake and weighted shake varieties; 

Fig. 10 shows stability levels for SHAKE-WT; 

Fig. 11 shows average delay for PIM in regular and shake varieties; 

Fig. 12 shows average delay for SLIP in regular, shake, and hybrid varieties; 

Fig. 13 shows average delay for RND in regular, shake, and hybrid varieties; 
40 Fig. 14 shows average delay for RR-RND in regular and shake varieties; 

Fig. 15 shows average delay for SHAKE; 

Fig. 16 shows average delay for ILPF in regular, shake, and weighted shake varieties; 

Fig. 17 shows average delay for RR-LQF regular, and weighted shake varieties; 

Fig. 18 shows average delay for SHAKE WT; 
45 Fig. 1 9 shows average delay for MBM heuristics with bursty traffic; 

Fig. 20 shows average delay for MWBM heuristics with bursty traffic; 

Fig. 21 shows average delay for MBM heuristics with different switch sizes; 

Fig. 22 shows average delay for MWBM heuristics with different switch sizes; 

Fig. 23 shows a flowchart depicting a preferred embodiment of the unweighted shakeup method; 
so Fig. 24 shows a flowchart depicting a preferred embodiment of the weighted shakeup method; 

Fig. 25 shows a preferred embodiment of an unweighted shakup switch; 

Fig. 26 shows a preferred embodiment of a weighted shakup switch; and 

Fig. 27 shows a preferred embodiment of a network system with weighted and unweighted shakeup switches. 

55 a. The Shakeup Technique 

[0046] The shakeup approach comes in two varieties: one for unweighted and another for weighted heuristics. The 
motivation behind the Shakeup operation in both cases is to provide a way for the scheduler to escape from a local 
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maximum, i.e. a maximal matching, and by doing so to speed the convergence of the queuing system to a steady state. 
Randomization is a crucial component of Shakeup. For the present discussions on bipartite graphs, let A be the set of 
vertices representing input ports, and let B be the set of vertices representing output ports. 

5 A.1 Unweighted Shakeup 

[0047] The preferred embodiment of the unweigted shakeup method is described using a flowchart shown in FIG. 
23 that describes the steps in the unweighted shakeup method. The switch used in this case has N input ports and N 
output ports. Each input port has N queues associated with it, each queue holding cells to be routed for a specific 

10 output port. Thus there are N x N queues. Cells to be transmitted are thus queued in the appropriate on of the N x N 
queues. At 23.1 each input port is represented as a first set of vertices in a bipartite graph. At 23.2 each output port is 
represented as a second set of vertices in the bipartite graph. At 23.3 non-empty queues are represented as edges in 
this bipartite graph. It should be clear that at a maximum there are N x N edges possible if ail the queues are non- 
empty. Some initial matching is then performed in 23.4; the initial matching has perhaps been generated through some 

15 conventional heuristic. In this initial matching a subset of edges is selected such that no vertex is incident to more than 
one edge. 

[0048] An unmatched vertex is now selected at 23.6. Note that the process loops for each edge in step 23.5. In a 
key step of the unweighted shakeup technique, the basic idea is that each unmatched vertex in A will be allowed to 
force a match for itself, even if that means removing an existing match. Specifically, each unmatched vertex x in set A 

20 that is incident to at least one edge selects uniformly at random one of those incident edges, say (x.y) in 23.7. A test 
is performed at 23.8 to see a match (x,y) exists. If y G B is unmatched, then (x, y) is added to the matching at 23.9. If 
y is already matched through edge (z, y), z e A, then (z, y) is removed from the matching and (x, y) is added at 23.11 . 
If more than one vertex in A attempts to match to the same vertex in B during a Shakeup, then one of the competing 
vertices in A is selected in a random fashion as the winner. Shakeup is not implemented in parallel but instead un- 

25 matched vertices are processed in A in a randomly selected order. Note that unweighted Shakeup returns a matching 
that is at least as large as the initial matching. 

[0049] FIG. 25 shows a preferred embodiment of an unweighted shakeup switch according to the present invention. 
This switch has four input ports 25.11-25.14 and four output ports 25.51-25.54. Cells are routed between the input 
ports and the output ports. An unweighted scheduler 25.6 performs the routing. The unweighted scheduler comprises 

30 a graph generator 25.2 that generates a bipartite graph with each input ports 25.11-25.14 as a node in a first set of 
vertices of the bipartite graph and each output ports 25.51-25.54 is a node in a second set of vertices of the bipartite 
graph. A specific cell to be routed between an input port and an output port is an edge between a vertex in the first set 
and the second set. The unweighted scheduler further comprises a matcher 25.3 for performing an initial matching by 
finding a subset of edges such that no vertex is incident to more than one edge in such subset. An unweighted shakeup 

35 matcher 25.4 performs an unweighted shakeup matching for a hitherto unmatched vertex x in said first set to a vertex 
y in the second set. The unweighted shakeup matcher adds new matching (x,y) if y is not part of any previous matching 
and replaces g a matching involving y by a matching (x,y) if y is already part of a matching (z,y). 

A.2 Weighted Shakeup 

40 

[0050] The preferred embodiment of the weigted shakeup method is described using a flowchart shown in FIG.24 
that describes the steps in the unweighted shakeup method. 

[0051] The switch used in this case has N input ports and N output ports. Each input port has N queues associated 
with it, each queue holding cells to be routed for a specific output port. Thus there are N x N queues. Cells to be 

*5 transmitted are thus queued in the appropriate on of the N x N queues. At 24.1 each input port is represented as a first 
set of vertices in a bipartite graph. At 24.2 each output port is represented as a second set of vertices in the bipartite 
graph. At 24.3 non-empty queues are represented as edges in this bipartite graph. It should be clear that at a maximum 
there are N x N edges possible if all the queues are non-empty. Some initial matching is then performed in 24.4; the 
initial matching has perhaps been generated through some conventional heuristic. Weighted Shakeup is slightly dif- 

50 ferent in that it takes into account queue lengths. The basic idea behind extending unweighted to weighted is to give 
preference to those input and output polls that are heavily loaded; the unweighted Shakeup applied above does not 
do this. Indeed, as all edges are treated equally in the unweighted Shakeup approach, it is possible that applying 
unweighted Shakeup to a weighted heuristic will make the solution worse. Weighted Shakeup, therefore, uses weights 
when selecting possible edges to add to the matching. Specifically, each unmatched vertex in set A that is incident to 

55 at least one edge randomly selects one of those edges in 24.6. Note that the process loops for each edge in step 24.5. 
This random selection of edges is done with a probability proportional to its relative size over ail the destination queues. 
That is, for input port i, if the queue holding cells destined for output j has length qy, then edge (i,j*) is selected with 
probability qj/(£k-i Nc lik)- Let (i,j') be the selected edge. If the vertex j' G B is unmatched, then (i,j') becomes part of the 
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matching 24.9. If j' is already matched through edge (kj 1 ), then with probability q^(qg, + q^,), (k,j') is removed from the 
matching and (i,j') is added 24.11 . Otherwise, the matching is not modified. 

[0052] Fig. 26 shows a preferred embodiment of an unweighted shakeup switch according to the present invention. 
This switch has four input ports 26.11-26.14 and four output ports 26.51-26.54. Cells are routed between the input 

5 ports and the output ports. A weighted scheduler 26.6 performs the routing. The weighted scheduler comprises a graph 
generator 26.2 that generates a bipartite graph with each input ports 26.11 -26.14 as a node in a first set of vertices of 
the bipartite graph and each output ports 26.51-26.54 is a node in a second set of vertices of the bipartite graph. A 
specific cell to be routed between an input port and an output port is an edge between a vertex in the first set and the 
second set. The weighted scheduler further comprises a matcher 26.3 for performing an initial matching by finding a 

10 subset of edges such that no vertex is incident to more than one edge in such subset. The weighted scheduler further 
comprises a weighted shakeup matcher 26.4 for performing a weighted shakeup matching for a hitherto unmatched 
vertex x in said first set to a vertex y in the second set, wherein said weighted shakeup matcher adds (x,y) if y is not 
part of any previous matching, said vertex y being selected with a probability: 

15 

gxy 
N 



wherein is a length of a queue between vertices x and y, and replaces a matching involving y by a matching (x,y) 
if y is already part of a matching (z,y), said replacement being done with a probability: 



[0053] Shakeups can be used in a stand-alone manner, or in conjunction with other heuristic approaches. Two general 
30 types of implementations are envisioned: orthogonal and hybrid. For the orthogonal type, the "basic heuristic* (e.g., 
SLIP, PIM, etc.) and Shakeup are performed in alternate phases. For example, if SLIP is the basic heuristic sequence 
of operations would be: SLIP, Shakeup, SLIP, Shakeup, SLIP Shakeup, etc. For the hybrid type, the Shakeup phase 
is performed simultaneously with the basic heuristic phase. Each unmatched input port examines which output ports 
it can send cells to. If any of these output ports is not matched, the input port competes for it as defined by the basic 
35 heuristic, and does not do Shakeup in this iteration. On the other hand, if all the desired output ports are already 
matched, then the input port performs Shakeup this iteration instead of the basic heuristic. Although the orthogonal 
and hybrid approaches are slightly different, the basic idea of using randomization to get out of local maxima is the 
same, and in practice both approaches are likely to have similar performance. For some implementations, however, 
the hybrid approach can require less communication than the orthogonal approach. In terms of physical implementation, 
40 the Shakeup schemes would be fairly simple to implement, using hardware similar to that needed by many of the 
existing unweighted heuristics. Some additional logic for a pseudo-random generator would be necessary, but it is 
estimated that this logic does not have to be sophisticated. The present simulations obtained good performance by 
implementing Shakeup using the C built-in pseudo-random generator. 

45 A 3 Network System with Shakeup Switches and Unshakeup Switches 

[0054] Fig. 27 shows a preferred embodiment of a network system with shakeup and unshakeup switches. This 
network system comprises a network 27.1 that comprises several switches. Of the switches, 27.31 and 27.32 are 
unshakeup switches and 27.20, 27.21 and 27.22 are unshakeup switches. 

so 

A.4 Theoretical Justification 

[0055] Theoretical justification for the present approach is discussed. This can be found in earlier work on bipartite 
matching. 

55 [0056] At first sight one may argue that there is no apparent reason why unweighted Shakeup should improve sched- 
uler performance and even converge to a maximal matching. Indeed, an unmatched input i sends only one request to 
an output j' at each iteration, and j' could very well be matched, in which case the size of the matching will remain the 
same. However by knocking off the matching the mate of j', a chain of knockoffs is triggered which leads, in a probabilistic 
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sense, to the discovery of an augmenting path-a standard technique for increasing the size of a given non-maximum 
matching. See C. Berge, "Two theorems in graph theory' Proceedings of the National Academy of Science, vol. 43, 
pp.842-844, 1957. This idea has been used with success in the design of probabilistic methods. Some basic facts are 
now provided. 

s [0057] Formally, a random bipartite graph G p (2N) is a graph with N vertices on each side of the bipartition, such that 
an edge between any two vertices on different sides is included with probability p independently of the others. The 
random graph model has been extensively studied, in part to capture the behavior of graph algorithms on lypicaP 
inputs. See B. Bollobas, Random Graphs. Harcourt Brace Janovich, 1 985 for an excellent survey of the subject. In the 
present problem p could be related to, although it is clearly different from, the probability that input i has a cell to send 

10 to output j for any pair (i,j) at the given time step at which a permutation needs to be found. It should be clear though 
that the bipartite graph of the switch cannot be accurately modeled by G p (2N), even for uniform traffic, due to the 
existence of queues. Consider a simplified unweighted Shakeup where in every iteration the unmatched vertices in A 
are considered in a fixed order. It was shown by Angluin and Valiant that this simplified Shakeup delivers with high 
probability a perfect matching on G p (2N) in 0(N log N) time, when p > (clog N)/N, for some constant c> 0. See D. 

is Angluin and L. C. valiant, "Fast probabilistic algorithms for Hamiltonian circuits and matchings," Journal of Computer 
and System Sciences, vol. 18, pp.155-193, 1979. Therefore, repeated application of unweighted Shakeup will give 
optimal throughput at a single time slot, under appropriate probabilistic assumptions on the traffic. This theoretical 
evidence is complemented by the experimental results in Section C. 

20 B. Framework for Comparison 

[0058] The use of the approach of the present invention is demonstrated herein with a variety of MBM and MWBM 
heuristics. As can be seen, in several cases, the range of switch stability increases from approximately 85% to over 
95% when the techniques of the present invention are used. It can also be seen that even for stable flows, the use of 

25 Shakeup can greatly reduce the average cell delay. 

[0059] Through extensive simulation, a number of different scheduling techniques are compared. 
[0060] It is believed that he range of stability is a fundamentally important metric for scheduling methods. An exper- 
imental definition of stability is provided herein and stability tests are run for a wide range of scheduling approaches, 
traffic types, and traffic loads. A number of different traffic types, designed to demonstrate the strengths and weak- 

30 nesses of scheduling techniques, are compared. The establishment of a benchmark suite of traffic types is an important 
contribution for the analysis of scheduling methods, particularly when the scheduling methods or traffic models are 
mathematically intractable. In particular, "low-degree" traffic flows, which substantially degrade the performance of 
MBM- and MWBM-related heuristics is compared. 

[0061] To compare the schedulers a variety of traffic types are also compared. In general, the goal is to begin to 
35 define a benchmark suite of traffic types that could be used to compare different schedulers. These traffic types should 
be selected to represent both typical and worst-case situations. 

[0062] For these simulations, each input port receives cells generated from its own source (that is, there are N 

sources). Each source is independent of the other sources, and generates at most one cell per time step. 

[0063] A variety of probability density functions are used for selecting output port destinations. One of the major 

40 contributions of this paper is the present investigation of low-degree traffic patterns, which incur poor performance for 
many schedulers. In this context, low-degree implies that the input and output ports of the bipartite graph have low 
degree. Intuitively, greedy heuristics should perform better for high-degree graphs than for low-degree graphs. High- 
degree graphs will have many good (i.e., optimal or near-optimal) matchings to choose from, while far fewer good 
matchings will exist for a low-degree graph. For a high-degree graph, a wealth of good solutions exist that greedily 

45 adding an edge to a matching probably will not greatly affect the solution quality. This is not the case for a low-degree 
graph, where early bad decisions can easily hurt the solution quality. 
[0064] Three types of destination probability density functions are also considered. 

B.1 High-degree balanced 

50 

[0065] Output ports are chosen with a uniform distribution over the N possibilities. All input ports and output ports 
will handle the same traffic load. 

B.2 Low-degree balanced 

55 

[0066] Each input port is the source of exactly k flows and each output port is the destination of exactly k flows. All 
flows have the same rate. To create the flows at the start of the simulation, k random permutations are used, selecting 
each of the N! possibilities with equal probability. 
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[0067] For example, let N = 8 and k = 2. Suppose the first permutation is (1 ,3,7,5,6,4,2,0) and the second permutation 
is (5,6,2,3,1,4,0,7). This implies input port 0 will send cells only to output ports 1 and 5 (with equal probability); input 
port 1 will send cells only to output ports 3 and 6; etc. Note that it is possible for an input port to have less than k 
possible destinations. In this example, input port 5 sends only to output port 4. 
s [0068] For this type of traffic, all input ports and output ports will have the same load. Also note that for the experiments 
described in Section C, only a single random set of permutations is used. 

B. 3 Low-degree unbalanced 

10 [0069] McKeown et al provide an example of an admissible flow pattern that causes the queues to grow unboundedly 
if MBM is the scheduling algorithm. See N. McKeown, V Anantharam, and I. Walrand, "Achieving 100% throughput in 
an input-queued switch," in Proceedings IEEE INFOCOM'96, (San Francisco, CA), pp. 296-302, March 1996. The 
example is essentially a case of unbalanced flows. This class is interesting because, intuitively, it is better to give 
preference to heavily loaded ports. For the present simulations, half of the input ports are the source of exactly k flows 

is and half are the source of 2k flows. Each flow generates traffic at the same rate, so half the input ports are loaded at 
rate X while the other half are loaded at 2 X. To create the flows at the start of the simulation, 2k random permutations 
are used, selecting each of the N! possibilities with equal probability. For the heavily loaded input ports, alt 2k flows 
are defined by the permutations. For the lightly loaded ports, only the first k permutations are used to define the flows. 
Note that the output ports will support no less than k flows and no more than 2k flows. 

20 [0070] A rate of x% implies all input and output ports are loaded at that rate. For low-degree unbalanced traffic, a 
rate of x% implies the heavily loaded input ports run at x%, while the lightly loaded input ports run at rate xJ2 %. 
[0071] The burstiness of the traffic can also effect performance. These two types are now considered: 
[0072] Non-bursty (Bernoulli): At each time step, a cell is generated by a source with some fixed probability, inde- 
pendent of all previous traffic. The probability density function that determined the output port is the same for all cells 

25 generated by the source. 

[0073] Bursty (Geometric): Traffic is generated as bursts of cells (ON bursts), followed by bursts of no cells (OFF 
bursts). The probability density function that determines the output port is the same for all ON bursts generated by the 
source. All cells within an ON burst will have the same destination. To create a burst with expected length m, the burst 
is discontinued at each time step with probability 1/m. (The use of a geometric distribution to model bursty traffic can 

30 be found in See N. W. McKeown, Scheduling Algorithms for Input-Queued Cell Switches. PhD thesis, University of 
California at Berkeley, 1995, for example.) 

[0074] Most of the present experiments assume non-bursty traffic. It will be seen see later that the experiments 
confirm empirical findings by McKeown that suggest a linear relationship between expected burst length and expected 
cell delay for a range of scheduling approaches See N. W. McKeown, Scheduling Algorithms for Input-Queued Cell 
35 Switches. PhD thesis, University of California at Berkeley, 1995. 

C. Experimental Results 

[0075] The experiments provided some valuable insights into the relative performance of different methods. In Table 
40 1 provide a summary of the experimental results are provided for the three most powerful methods MBM, MWBM-LQF 
and MWBM-LPF. The performance of these methods should be the benchmark against which heuristics are tested. It 
is interesting to notice how close the performance of MBM is to the performance of the weighted matching methods 
for the inputs that are examined. Ideally, unweighted heuristics should perform well also, and often they do, but it is 
demonstrated that there is considerable room for improvement. 

45 

Table 1 



Scheduler 


Phb 


P tb 


Plu 


Phb 


Pib 


P lu 


MBM 


100% 


100% 


97% 


2.90 


2.86 


1.47 


MWBM- 


100% 


100% 


100% 


2.94 


2.70 


1.35 


LQF 














MWBM- 


100% 


100% 


100% 


2.87 


2.23 


1.15 


LPF 














Table 1 : Experimental stability rates (p) and average delays (5) for MBM and MWBM approaches. High-degree 
balanced (hb), low-degree balanced (lb), and low-degree unbalanced (lu) flows are used for 32 ports. Average 
delays are for 100,000 time steps with 85% traffic load. 
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[0076] A general conclusion that can be drawn from the test data is that the addition of Shakeup or weighted Shakeup 
on top of most of the examined heuristics improves performance. Another lesson from the experiments is that most of 
the previously known heuristics reach their full potential in terms of stability typically in 2-3 iterations. However Shakeup 
typically exhibited slowly improving performance as the number of iterations increased beyond this point. The exper- 
s iments were designed to fall into three major categories. Each category is elaborated herein providing the description 
of the different settings and specific comments on the observed results. 

C.1 Stability vs. Iterations 

10 [0077] In the first category the various methods are tested to determine the maximum input rate for which they are 
stable as a function of the number of iterations of the method to select a single permutation. This is called the stability 
rate. 

[0078] The present definition of stability is by necessity experimental. Let A be a method running for i iterations. The 
range of traffic rates from 0% to 100% is considered, as defined for the traffic types described in Section B. For every 

is rate value r queried by the search mechanism, the system is simulated for at most T = 256N 2 steps. At every time step, 
a check is done to see if the sum of all the queue lengths Q has become greater than or equal to 1 6N 2 if so, the method 
is declared unstable at rate r. Also, checkpoints are established for every 1,000 time steps; if the value of 0 at a 
checkpoint is equal to or less than the value of Q at the previous checkpoint, the method is declared stable at rate r. 
Finally, if the method is not declared either stable or unstable after T time steps, the method is declared stable at rate 

20 r . If the method is stable for rate r, the search mechanism considers only larger rates in the future. If the method is 
unstable, only smaller rates are considered in the future. For efficiency, in the actual implementation a linear search 
is not performed on all possible rates but a binary search is performed instead. 

[0079] In summary, the experimental space for the first category of experiments consists of all the possible values 
of the tuple (A, i, r, flow type). The number of iterations i ranges from 1 to 8. Observe that for some of the approaches 

25 (MEM, RR-RND, ILPE and RR-LQF), iterations are meaningless since maximal matchings are computed in the first 
iteration. For those methods a horizontal line is depicted in the graph. The two axes show rate and number of iterations. 
The relevant graphs are in Figures 3-10. In each graph for this category the results for related methods A are aggre- 
gated. Typically there will be three methods in a graph, called X, X-HYBRID and X-SHAKE. X is any vintage method 
like PIM or SLIP. X-SHAKE, implements the orthogonal use of Shakeup: at every one of the i iterations X is run followed 

30 by an execution of unweighted Shakeup, X-HYBRID, implements the hybrid use of Shakeup: at each of the i iterations 
the basic heuristic X 'merged' with Shakeup as outlined in Section A is run. 

C.1 (a) Unweighted methods 

35 [0080] The graphs for this group of experiments demonstrate that the addition of Shakeup to the unweighted heu- 
ristics can typically improve the maximum stability rate by 10%. The increase becomes palpable after already two 
iterations and slowly increases thereafter with the number of iterations. As expected, the increase is smaller for the 
"easy" high-degree balanced input, on which almost all unweighted heuristics achieve by themselves stability rate of 
95% or higher after 3 iterations. 

40 

C.1 (b) Weighted methods 

[0081] The graphs for the I LPF method (Fig. 8) show that Shakeup provides more modest improvements than in the 
unweighted case, of the order of 2-3%. This is not surprising given than ILRF is a heavy-duty centralized method that 
45 does global sorting. 

[0082] As expected, weighted Shakeup (SHAKE-WT) performs slightly better than pure Shakeup which can margin- 
ally degrade performance every other iteration (due to the fact that it blindly unmatches the heavily loaded ports matched 
by ILPF). The situation is similar for RR-LQF although the gains from (weighted) Shakeup are larger, of the order of 
5-6%. 

50 [0083] Finally, to demonstrate the robustness of Shakeup and weighted Shakeup are tested as stand-alone methods. 

The results are shown in Figures 7 and 1 0. The results are particularly good given the very simple nature of the methods. 

If enough (i.e. 8) iterations are performed, both methods achieve stability rate of 90% or higher on all types of inputs. 

In the low-degree unbalanced input, which was specifically designed to be harder for unweighted methods, weighted 

Shakeup achieves a stability rate of 98% after 6 iterations while Shakeup is barely above 90% at the same number of 
55 iterations. 
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C.2 Average Delay vs. Iterations 

[0084] Average delay means the average delay of all cells that were successfully routed within the time window of 
the simulation. (In general, the input queues can still contain cells that do not contribute to this value.) A cell that is 

s routed immediately has delay 0. The stability experiments outlined above, give information that allows in most cases 
a clear ranking of the methods based on the maximum input rate they can support. However, this ranking would not 
give the full picture. For a given rate r it could very well be the case that both methods A and A' are declared (eventually) 
stable but they differ significantly in the average delay they incur on the packets. On the other hand, testing the methods 
solely based on average delay and ranking them accordingly could be misleading as well. First, a method may actually 

io settle for a large but stable average delay for a given rate r, information that is lost by looking at the average delay 
only. Second, what is the rate r at which one should examine the average delay incurred by A? If r is high enough so 
that A is unstable, the average delay obtained will be essentially meaningless and A will be judged unfairly under inputs 
it cannot possibly handle. Using the data from the stability results this pitfall is avoided by choosing r = 85% for all the 
average delay graphs, a value at which almost all of the methods tested were shown to be stable. In summary, the 

is experimental space for the Avenge Delay Graphs consists of the possible values of the tuple (A, i, r = 85%, flow type). 
The two axes show avenge delay (in time steps) and number of iterations. The relevant graphs are in Figures 11 -18. 
[0085] The average delay graphs confirm the picture given by the stability graphs for the advantages gained by 
Shakeup. The gains are more tangible for the unweighted methods. The improvements obtained by the application of 
Shakeup range from roughly 10% to a quite dramatic 10-fold improvement for SLIP and RR-RND on low-degree bal- 

20 anced inputs. Again in Figures 1 5 and 1 8 Shakeup and weighted Shakeup are tested as stand-alone methods to verify 
robustness. After 6 iterations, the average delay is competitive with the one obtained by almost any other heuristic. 

C.3 Average Delay vs. Burst Size or Ports 

25 [0086] In this final category, the qualitative exploration of the different methods are continued by plotting average 
delay against two other parameters of importance. 

[0087] In the first case, the experimental space consists of all possible values of the tuple (A, i = 9, r = 85%, burst 
size, flow type). The axes show average delay and burst size. See Figures 1 9 and 20. The number of time steps is set 
again at 100,000. The burst size increases by increments of 17 in the range [See A. C. Kam and K.-Y, Siu, "Linear 

30 complexity algorithms for bandwidth reservations and delay guarantees in input-queued switches with no speedup," 
in Proceedings of the 6th International Conference on Network Protocols, pp. 2-1 1 , 1 998 and 1 36]. The results are as 
expected, i.e. for all methods the average delay increases as a function of the burst size, and this increase happens 
essentially at the same rate for all the methods. The curve for the heavy-duty MBM and MWBM methods is essentially 
the same as the ones for the rest, which implies that the increase of the average delay with the burst size is essentially 

35 unavoidable. The conclusion is that increased buffer space is necessary to store packets if bursty traffic is anticipated. 
[0088] In the second case, the experimental space consists of all possible values of the tuple (A, s = 9, r = 85%, p, 
flow type) where p denotes number of ports. The axes show average delay and number of ports. For p = 4, 8,16, 32, 
64. Figures 21 and 22 are tested and the results for unweighted- and weighted-type methods are presented respectively. 
The results demonstrate that most of the methods scale reasonably well, i.e. the average delay increases very slowly 

40 as the number of ports grows. There are some exceptions though. The average delay for RR-RND in Fig. 21 increases 
sharply with the number of ports for low-degree balanced flows. RR-RND is by definition a very rigid method in the 
way it prioritizes input ports. As the number of ports increases, this inflexibility incurs a higher overhead. The average 
delay for SLIP in Figure 21 fluctuates between 30 and 140 for low-degree balanced and between 3 and 26 for unbal- 
anced inputs. SLIP'S behavior is hard to predict a priori due to the different priority wheels maintained by each port. If 

45 these wheels are aligned the wrong way - an event which apparently occurs on the same traffic type simply as N grows 
- SLIP'S performance suffers. The variation in average delay for the weighted-type methods in Fig. 22 is bounded by 
5; therefore these methods scale in a satisfactory manner. 

[0089] Other modifications and variations to the invention will be apparent to those skilled in the ad from the foregoing 
disclosure and teachings. Thus, while only certain embodiments of the invention have been specifically described 
so herein, it will be apparent that numerous modifications may be made thereto without departing from the spirit and scope 
of the invention. 



Claims 

55 

1 . A computer-implemented unweighted queuing method for a system with M groups of N input queues each of said 
N input queues requesting a unique one of N services, wherein M and N are natural numbers, characterized by 
the steps of: 
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a) (23.1) representing each of said M group as a node in a first set of vertices of a bipadite graph; 

b) (23.2) representing each of said N services as a node in a second set of vertices of said bipartite graph; 

c) (23.3) representing each non-empty queue as an edge between a vertex in said first set corresponding to 
one of said M groups to which said each non-empty queue belongs to and a vertex in said second set corre- 
sponding to a service requested by said each non-empty queue; 

d) (23.4) performing an initial matching by finding a subset of edges such that no vertex is incident to more 
than one edge in such subset; 

e) (23.6-23. 1 1 ) performing an unweighted shakeup matching for a hitherto unmatched vertex x in said first set 
to a vertex y in the second set by performing an appropriate one of following two steps, wherein if x matches 
with more than one of said vertices in the second set, vertex y being selected at random: 

e-1)(23.9) adding new matching (x,y) if y is not part of any previous matching; and 

e-2)(23.11) replacing a matching involving y by a matching (x,y) if y is already part of a matching (z.y), 

wherein z is a vertex in said first set; and 

f) (23.5, 23.10) repeating step (e) for each unmatched vertex in said first set. 

A computer-implemented weighted queuing method for a system with M groups of N input queues each of said N 
input queues requesting a unique one of N services, wherein M and N are natural numbers, characterized by the 
steps of: 

a) (24.1) representing each of said M groups as a node in a first set of vertices of a bipartite graph; 

b) (24.2) representing each of said N services as a node in a second set of vertices of said bipartite graph; 

c) (24.3) representing each non-empty queue as an edge between a vertex in said first set corresponding to 
one of said M groups to which said each non-empty queue belongs to and a vertex in said second set corre- 
sponding to a service requested by said each non-empty queue; 

d) (24.4) performing an initial matching by finding a subset of edges such that no vertex is incident to more 
than one edge in such subset, wherein if a vertex can be matched using more than one edges, one edge is 
selected from said more than one edges with a probability relative to the size of queues associated with each 
of said one or more edges; 

e) (24.6-24. 11) performing a weighted shakeup matching by performing an appropriate one of following two 
steps for a hitherto unmatched vertex x in said first set to a vertex y in the second set: 

e-1)(24.7) if y is not part of any previous matching, said vertex y being selected with a probability: 

Qxy 
N 

k=l 

wherein q^ is a length of a queue between vertices x and Y; 

e-2)(24.11) replacing a matching involving y by a matching (x,y) if y is already part of a matching (z.y), 
said replacement being done with a probability: 

wherein z is a vertex in said first set; and 

f) (24.5, 24.10) repeating step (e) for each unmatched vertex in said first set. 

A method of unweighted routing in a switch comprising input ports and output ports, having cells to route between 
said input ports and said output ports, each input port having one queue each for each said output port, thereby 
each input port having as many queues as a number of said output ports, characterized by the steps of: 
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a) representing each said input pods as a node in a first set of vertices of a bipartite graph; 

b) representing each said output ports as a node in a second set of vertices of said bipartite graph; 

c) representing each non-empty queue between any of said input ports and any of said output ports as an 
edge between a vertex in said first set corresponding to said any of said input ports and a vertex in said second 
set corresponding to said any of said output ports, thereby creating a set of edges; 

d) performing an initial matching by finding a subset of edges such that no vertex is incident to more than one 
edge in said subset; 

e) performing an unweighted shakeup matching for a hitherto unmatched vertex x in said first set to a vertex 
y in the second set by performing an appropriate one of the following two steps, wherein if x matches with 
more than one of said vertices in the second set, vertex y being selected at random: 

e-1 ) adding new matching (x,y) if y is not part of any previous matching; 

e-2) replacing a matching involving y by a matching (x,y) if y is already part of a matching (z,y) wherein z 
is a vertex in said first set; and 

f) repeating step (e) for each unmatched vertex in said first set. 

A method of weighted routing in a switch comprising input ports and output ports, having cells to route between 
said input ports and said output ports, each input port having one queue each for each said output port, thereby 
each input port having as many queues as a number of said output ports, characterized by the steps of: 

a) representing each said input ports as a node in a first set of vertices of a bipartite graph; 

b ) representing each said output ports as a node in a second set of vertices of said bipartite graph; 

c) representing each non-empty queue between any of said input ports and any of said output ports as an 
edge between a vertex in said first set corresponding to said any of said input ports and a vertex in said second 
set corresponding to said any of said output ports, thereby creating a set of edges; 

d) performing an initial matching by finding a subset of edges such that no vertex is incident to more than one 
edge in such subset, wherein if a vertex can be matched using more than one edges, one edge is selected 
from said more than one edges with a probability relative to a size of a queue associated with each of said 
one or more edges; 

e) performing a weighted shakeup matching by performing an appropriate one of following two steps for a 
hitherto unmatched vertex x in said first set to a vertex y in the second set: 

e-1 ) if y is not part of any previous matching, said vertex y being selected with a probability: 



N 

wherein q^ is a length of a queue between vertices x and Y; 

e-2) replacing a matching involving y by a matching (x,y) if y is already part of a matching (z,y), said 
replacement being done with a probability: 

wherein z is a vertex in said first set; and 
f ) repeating step (e) for each unmatched vertex in said first set. 
The method of claim 3 wherein said unweighted shakeup is performed for every alternate matching. 
The method of claim 4 wherein said weighted shakeup is performed for every alternate matching. 
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7. The method of claim 3 wherein said unweighted shakeup is performed only if no increase in the size of matching 
is possible without shakeup. 

8. The method of claim 4 wherein said weighted shakeup is performed only if no increase in the size of matching is 
possible without shakeup. 

9. An unweighted shakeup switch for use in a network, said switch comprising: 

input ports (25.11-25.14); 

output ports(25.51 -25.54), wherein cells are routed between said input ports and said output ports; and 
an unweighted scheduler, 

wherein said unweighted scheduler further comprises: 

a graph generator (25.2) that generates a bipartite graph with each said input ports as a node in a first set of 
vertices of the bipartite graph and each said output ports is a node in a second set of vertices of said bipartite 
graph and a cell to be routed between an input port and an output port is an edge between a vertex in said 
first set said second set; 

a matcher (25.3) for performing an initial matching by finding a subset of edges such that no vertex is incident 
to more than one edge in such subset; 

an unweighted shakeup matcher (25.4) for performing an unweighted shakeup matching for a hitherto un- 
matched vertex x in said first set to a vertex y in the second set, wherein said unweighted shakeup matcher 
adds new matching (x,y) if y is not part of any previous matching and replaces a matching involving y by a 
matching (x,y) if y is already part of a matching (z,y). 

10. A weighted shakeup switch for use in a network, said switch comprising: 

input ports (26.11-26.14); 

output ports (26.51-26.54), wherein cells are routed between said input ports and said output ports; and 
a weighted scheduler, 

wherein said weighted scheduler further comprises: 

a graph generator (26.2) that generates a bipartite graph with each said input ports as a node in a first set of 
vertices of the bipartite graph and each said output ports is a node in a second set of vertices of said bipartite 
graph and a cell to be routed between an input port and an output port is an edge between a vertex in said 
first set said second set; 

a matcher (26.3) for performing an initial matching by finding a subset of edges such that no vertex is incident 
to more than one edge in such subset; 

a weighted shakeup matcher (26.4) for performing a weighted shakeup matching for a hitherto unmatched 
vertex x in said first set to a vertex y in the second set, wherein said weighted shakeup matcher adds (x,y) if 
y is not part of any previous matching, said vertex y being selected with a probability: 

Qxy 
N 

wherein is a length of a queue between vertices x and y, and replaces a matching involving y by a matching 
(x,y) if y is already part of a matching (z,y), said replacement being done with a probability: 

11. The switch of claim 9 wherein said unweighted shakeup is performed for every alternate matching. 
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12. The switch of claim 10 wherein said weighted shakeup is performed for every alternate matching. 

13. The switch of claim 9 wherein said unweighted shakeup is performed only if no increase in the size of matching 
is possible without shakeup. 

s 

14. The switch of claim 10 wherein said weighted shakeup is performed only if no increase in the size of matching is 
possible without shakeup. 

15. A network system comprising a plurality of switches for routing messages, a first subset of the plurality of said 
10 switches being unweighted shakeup switches, a second subset of the plurality of said switches being weighted 

shakeup switches, 

each of the unweighted shakeup switches comprising: input ports (25.11-25.14); output ports(25.51 -25.54), 
wherein cells are routed between said input ports and said output ports; and an unweighted scheduler, wherein 
is said unweighted scheduler further comprises: 

a graph generator (25.2) that generates a bipartite graph with each said input ports as a node in a first set of 
vertices of the bipartite graph and each said output ports is a node in a second set of vertices of said bipartite 
graph and a cell to be routed between an input port and an output port is an edge between a vertex in said 
first set said second set; 

20 a matcher (25.3) for performing an initial matching by finding a subset of edges such that no vertex is incident 

to more than one edge in such subset; and 

an unweighted shakeup matcher (25.4) for performing an unweighted shakeup matching for a hitherto un- 
matched vertex x in said first set to a vertex y in the second set, wherein said unweighted shakeup matcher 
adds new matching (x,y) if y is not part of any previous matching and replaces a matching involving y by a 
25 matching (x,y) if y is already part of a matching (z,y), and 

each of the weighted shakeup switches comprising: input ports (26.11-26.14); output ports (26.51-26.54), 
wherein cells are routed between said input ports and said output ports; and a weighted scheduler, wherein 
said weighted scheduler further comprises: 

a graph generator (26.2) that generates a bipartite graph with each said input ports as a node in a first set of 
30 vertices of the bipartite graph and each said output ports is a node in a second set of vertices of said bipartite 

graph and a cell to be routed between an input port and an output port is an edge between a vertex in said 
first set said second set; 

a matcher (26.3) for performing an initial matching by finding a subset of edges such that no vertex is incident 
to more than one edge in such subset; and 
35 a weighted shakeup matcher (26.4) for performing a weighted shakeup matching for a hitherto unmatched 

vertex x in said first set to a vertex y in the second set, wherein said weighted shakeup matcher adds (x,y) if 
y is not part of any previous matching, said vertex y being selected with a probability: 



40 



45 



SO 



fry 
N 

Jt=l 

wherein is a length of a queue between vertices x and y, and replaces a matching involving y by a matching 
(x,y) if y is already part of a matching (z,y), said replacement being done with a probability: 
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FIG.6 
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FIG. 7 
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FIG.8 
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FIG. 9 
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FIG. 10 
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FIG. 12 
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FIG. 13 



AVERAGE DELAY, HIGH-OEGREE BALANCED DISTRIBUTION, 
32 PORTS. 85% RATE. 100000 STEPS 




RND-HYBRID--+- 
RNO-SHAKE--Q- 



^-,J=4^£^ &- 



4 5 6 7 B 

ITERATIONS 

AVERAGE DELAY, LOW-DEGREE BALANCED DISTRIBUTION. 
M PORTS. 85% BATE. 100000 STEPS 




BHD- 
RND-HYBRID- 
RNO-SHAKE- 



Vi T 

D -o 



• • • • -f^V • ~~~~ - - ~ ARt |4j*-*i> -a 



f 2 3 4 5 6 7 8 

ITERATIONS 

AVERAGE DELAY, LOW-DEGREE UNBALANCED DISTRIBUTION, 

32 PORTS, 85% RATE. 100000 STEPS 
100008 — — — — 




RND-+- 
RND-HYBRID--+- 
RND-SHAKE-O- 



— 4 — 



3 4 5 6 

ITERATIONS 



30 



EP 1 049 352 A2 



FIG. 14 

AVERAGE DELAY. HIGH-DEGREE BALANCED DISTRIBUTION. 
32 PORTS. 85% RATE. 100000 STEPS 
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FIG. 16 



AVERAGE DELAY, HIGH-OEGREE BALANCED DISTRIBUTION. 
32 PORTS. 85% RATE. 100000 STEPS 
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AVERAGE DELAY, HIGH-DEGREE BALANCED DISTRIBUTION, 
32 PORTS. 85% RATE. 100000 STEPS 
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AVERAGE DELAY. HIGH-DEGREE BALANCED DISTRIBUTION, 
32 PORTS. 85% RATE. 100000 STEPS 
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AVERAGE DELAY, HIGH-DEGREE BALANCED DISTRIBUTION. 
32 PORTS. 85% RATE. 100000 STEPS 
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